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The correlation function of the quantum fluctuations due 
to collective excitations is calculated and used to investigate 
the phase diffusion of a Bose-Einstein condensate close to zero 
temperature. It is shown that the phase diffusion time of the 
condensate is much longer than the result obtained by assum- 
ing that the correlation time of the quantum fluctuations is 
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I. INTRODUCTION 

The development of the technologies of laser trap- 
ping and evaporative cooling has yielded intriguing Bose- 
Einstein condensates (BECs) a state of matter in 
which many atoms are in the same quantum mechani- 
cal state. The remarkable observations of gaseous BECs 
have opened up new avenues of research into the 
physical properties and nature of Bose-condensed sys- 
tems. The phase properties of a BEC are of particular 
interest because the phase of an order parameter, i.e., the 
macroscopic wave function of the condensate, reflects di- 
rectly the coherent nature of the condensate. 

For BEC created in experiment, one of the most im- 
portant characters is that all the atoms in the conden- 
sate can be described by the wave function (i.e., the or- 
der parameter) with a single phase. Due to thermal and 
quantum fluctuations, however, the single phase of the 
condensate will become unpredictable beyond the phase 
diffusion time. After the realization of BECs, the phase 
diffusion of the condensate has been discussed intensively 
|Q ^0j. In particular, the role of quantum fluctuations on 
the phase diffusion process was investigated in the pio- 
neering work by Lewenstein and You J9|,[L0||, and a far 
off-resonant light scattering experiment was proposed to 
detect the quantum diffusion. Recently, a Langevin equa- 
tion was given by Graham [ ^4|Jl5|l to discuss the phase 
diffusion due to quantum fluctuations and thermal fluc- 
tuations. The calculation of the time scale of the phase 
diffusion is a very important problem because the phase 
diffusion time determines when the phase of a BEC would 
be unpredictable. Recently, the phase correlation has 
been investigated experimentally by the JILA group [^l) . 
It was found that there is no detectable diffusion of the 
phase on time scale 100 ms. The stable interference pat- 



terns shown in the experiment put forward a question |2l| ] 
why the phase correlation is so robust despite the phase 
diffusion and complicated rearrangement dynamics of the 
two condensates. 

In the present work, we address the question of the 
phase diffusion process of a condensate close to zero tem- 
perature. In general, the phase diffusion of the conden- 
sate can have either a thermal or a quantum origin. At 
extremely low temperature (in the experiment by the 
JILA group ipiJI ) the temperature is only 0.1T C , where 
T c denotes the critical temperature of the Bose gas.), the 
thermal fluctuations can be omitted and hence the quan- 
tum fluctuations become dominant. We give therefore 
emphasis on the role of collective excitations due to quan- 
tum fluctuations in a phase diffusion process. Although 
the phase diffusion process due to quantum fluctuations 
has been investigated by several authors such as the re- 
cent researches in |l4||l5|,^0| , the analysis of the time cor- 
relation of the quantum fluctuations is not given when 
the phase fluctuations are calculated. Obviously, the cor- 
rect consideration of the time correlation of the quantum 
fluctuations would make more reliable prediction on the 
phase diffusion process. In particular, researches show 
that the phase diffusion time is much longer than the cor- 
relation time of the quantum fluctuations. In this case, 
our results show that the phase diffusion time calculated 
from the correlation function of the quantum fluctuations 
is much longer than that obtained in the previous theo- 
retical researches |l4| , |l5| , p0[ . 

The paper is organized as follows. In Section II, we 
investigate the phase fluctuations of the condensate due 
to quantum fluctuations. In Sec. Ill, the phase diffusion 
time is calculated for the condensate close to zero tem- 
perature, where the effect of thermal fluctuations can be 
omitted. Sec. IV contains a discussion and summary of 
our results. 



II. PHASE FLUCTUATIONS OF THE 
CONDENSATE DUE TO QUANTUM 
FLUCTUATIONS 

For temperature below the critical temperature T c , the 
condensate can be described very well by the following 
order parameter with a phase factor <f> (t) 



$(r,i) =$ (r)e" 



i<p(t) 



(1) 



f 



where the phase of the condensate has the form 

4>{t)= l x{N< h T)t/U, (2) 

and the time-independent real component $o (r) is deter- 
mined by the stationary Gross-Pitaevskii (GP) equation 



H(N ,T) $ (r) 



-^V 2 + y cxt (r)+ 5 $g(r) U (r) 
2m J 

(3) 



where V^xt ( r ) is an external harmonic potential, and 
g = Airh 2 a s /m is the coupling constant fixed by the 
s— wave scattering length a s . The chemical potential 
fj, (Nq, T) in the above equation is determined by the nor- 
malization condition for the density distribution no ( r ) 
of the condensate. With a Thomas- Fermi approximation 
j| , one gets easily the following expression for the chem- 
ical potential: 
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(4) 
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where Who = (^x^y^z) is the geometric average of os- 
cillator frequencies, and ai 10 = y^Ti/mujho is the harmonic 
oscillator length of the system. From Eqs. (§) and (|J), 
we see that the particle number fluctuations of the con- 
densate yield fluctuations in the chemical potential, and 
hence lead to the phase diffusion of the condensate. 

Assuming the mean ground state occupation number 
is (N ) , the average phase of the condensate is then given 
by 



0((JVo),t) = A*«JVo),T)t/R. 



(•5) 



The phase diffusion of the condensate can be described 
by considering the phase difference A<j> (i) — cf>(t) — 
(j>((N ) ,t). From Eqs. (§) and (||), it is straightforward 
to obtain a differential equation on A(f> (t): 



dA(f> (t) 
Jt 



— A qua (t) — 



g> M ((A ),T) 
d(N ) 



AN (t) /h, (6) 



where A qua (t) is determined by collective excitations 
due to quantum fluctuations. In the above expression, 
ANo (t) represents the fluctuations of the ground state 
occupation number around (No). A similar equation was 
derived and used by Graham 0,^5) to discuss the phase 
diffusion of the condensate. Note that AAo (t) can be 
either negative or positive numbers. For AAo (t) < 0, 
there are collective excitations created so that the atoms 
would loss in the condensate. Similarly, AA (t) > 
means the annihilation of collective excitations, and the 
ground state occupation number would increase in this 
case. In addition, AAo (t) should be time-dependent be- 
cause it originates from quantum fluctuations. Eq. <M) 



is our starting point to discuss the phase diffusion of the 
condensate. Because it is obtained from Eqs. (||) and (||), 
rather than directly from a time-dependent GP equation, 
we anticipate that Eq. (||) is still correct for longer time 
where the time-dependent GP equation may be no longer 
valid |§. 

From Eq. (||) , the phase fluctuations of the condensate 
are given by 

((A<t> (r)) 2 ) = f f <A qua (0 A qua (£')) (7) 
Jo Jo 

where (A qua (£) A qua (£')) is the correlation function of 
the quantum fluctuations. When obtaining Eq. (0), we 
have assumed that ((A0) 2 ) = at time t = 0. The cal- 
culation of (A qua (£) Aq U a (£')) plays & crucial role in in- 
vestigating the phase diffusion process close to zero tem- 
perature. 

For the Bose gas trapped in a harmonic potential, it 
is convenient to use the following decomposition of the 
particle field operator 



{j, (f) = $ (r*) + ^ (ui (f) at + v* (r) at 



(8) 



where $ (f) = (tp (r)y is the well known order parameter, 
and the index i labels the elementary excitations of the 
system. For the collective excitations discussed here, the 
energy of the collective mode indexed by nl is given by 
the dispersion law p3[ 



e n i = TiL0h o (2n + 2nl + 3n - 



A/2 



(9) 



As shown in |2jg5(] , the contributions to condensate fluc- 
tuations due to quantum fluctuations are dominated by 
these phonon-type collective excitations. For the collec- 
tive mode nl, one obtains the following leading behaviour 
for u n i (r) and v nl (r) & 



u n i (r) ^ -Vni (r) 



I gno (r) 
2e„i 



Xm (r) , 



(10) 



where Xnl (^0 is the velocity potential associated 
with the collective mode, and satisfies the condition 
J drXni (r) Xij {r)d 3 f — <5 n i,ij. In addition, the average 
occupation number of the atoms corresponding to the 
collective mode indexed by nl is given by [f27|,[28| 



(Nnl) = Kl + O/nl: 



(11) 



where f n \ = [exp (e n \/kBT) — lp . When a collective ex- 
citation with index nl is created from the condensate due 
to quantum fluctuations, its energy e n \ originates from 
the energy fluctuations AE of the condensate. Under this 
consideration, a time-energy uncertainty relation can be 
used to calculate the longevity r n i of the collective mode 



2 



nl. The longevity r n \ of the collective mode nl is there- 

1/2 

fore approximated as l/who (2n 2 + 2nl + 3n + 1) . For 
JILA experiment ^J, this means that the longevity of 
the collective mode is smaller than 10 ms, which is obvi- 
ously much smaller than the phase diffusion time. 

When all collective modes are considered, A qua (t) can 
be written as: 



X, 



qua {t > ~ d (N ) 



nl^O 



(12) 



where AA„i (t) reflects the changes of the ground state 
occupation number due to the creation and annihila- 
tion of the collective mode nl. Therefore, the magni- 
tude of AN n \(t) can be regarded as (N n \). Note that 
AiV„i (t) itself can be either positive or negative, and 
varies with time due to quantum fluctuations. In the 
case of AN n \(t) < 0, there are |AA n i(i)| atoms cre- 
ated from the condensate due to quantum fluctuations, 
while AiV n ] (f) > represents the annihilation of AA n i (t) 
atoms. 

From Eq. (fjj), (A qua (£) A qua (£')) can be written as: 



(*qua (£) *qua (£')) = Qquae-| ? - 5 'l /7 



(16) 



When the above exponential form is used, r qua should be 
regarded as the average correlation time of the collective 
excitations, and is determined by the following expres- 



qua 



JZodTT 2 (X qua (t)JT qua (t + r)) 
J^driX^^X^it + r)) ■ 



(17) 



In terms of Eqs. (Q) and (|l7|), the correlation time of 
the quantum fluctuations is given by r qua = \/2/u\ w . In 
addition, in Eq. (|l6|), the magnitude Q qua of the corre- 
lation function is given by 

Qqua = 4 i d V i( N o) > T) jd (A )] 2 (S 2 N qu ,) . (18) 
n 

In the above expression, (£ 2 A qua ) = X)ni^p(^ni) 2 can be 
taken as the particle number fluctuations p5[ of the con- 
densate due to the collective excitations. Using the for- 
mulas (P|)-(|ll|), after a straightforward (although rather 
complex) calculation, we obtain the result of (5 2 iV qua ): 



<*qua (0*qua O 



^[^(<A ),T)/d(A )] 2 x 



£ E (AA nl (0AA n ,,,(O). (13) 

Assuming that there is no correlation between different 
collective modes, i.e., (AN n \ (£) AA n q< (£')) = when 
nl n'V , one gets the following expression for the corre- 
lation function: 



(A qua (0 X qua (£')) = -2 W {{No) , T) Id (A )] 2 x 
n 



^(A^e-M'l^ 

nl^O 



(14) 



When obtaining the above result, we have used the fol- 
lowing relation 



(AA nl (£) AA nl (O) = <A n 



(15) 



In the above expression, r n \ is the longevity of the col- 
lective excitation nl. When |£ — £'| is much larger than 
r n i, the correlation between the collective excitations at 
times £ and £' can be omitted. Therefore, r n i can be 
approximated as the correlation time of the correlation 
function (AA nl (£) AA nl (£')>■ 

Because the contributions to the quantum fluctuations 
come mainly from the low- lying collective excitations, 
as a reasonable approximation, the correlation function 
(A qua (£) A qua (£')) can b e approximated as an exponen- 
tial form: 



(o" 2 A qua ) = 0.958 



f 0,3 

V Oho 
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14.174 {-^-] A 12 / 15 , 
ah, 



4/5 



(19) 



where A is the total number of atoms in the trap. The 
second term on the right hand side of the above equa- 
tion represents the fluctuations due to the effect of the 
quantum depletion which is given in |24j] . This term has 
a finite contribution to the condensate fluctuations when 
the temperature approaches zero. Thus, we anticipate 
that there is still phase diffusion in the case of zero tem- 
perature. 

We now turn to discussing the phase diffusion of the 
condensate due to quantum fluctuations. From Eqs. (Q) 
and (fl6|), the phase fluctuations of the condensate, which 
play a crucial role in discussing the phase diffusion, read 



((A0(r)) 2 ) =2Q 



qua ' qua 



' qua 



' qua c 



,-t/t 



(20) 

The phase diffusion time r p hase can be obtained by setting 
{(A(f) (t)) 2 ) = 7r 2 in the above expression. 



III. PHASE DIFFUSION TIME OF THE 
CONDENSATE CLOSE TO ZERO 
TEMPERATURE 

We now turn to discussing the phase diffusion time 
using the phase fluctuations given by Eq. (p0|). It is 
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useful to discuss the phase fluctuations given by Eq. (J2C 
for two special cases. When the time r is much larger 
than the time scale of the correlation time r qua , the phase 
fluctuations of the condensate can be approximated as: 



((A</>(r)) 2 )«2Q 



qua ' qua ' 



(21) 



Therefore, if the phase diffusion time Tphaso calculated 
from Eq. ( |20| ) is much larger than r qua , the phase diffu- 
sion time in this situation takes the following analytical 
form: 



Tphasc 



2Q C 



(22) 



In the case of r << r qua , however, the phase fluctua- 
tions of the condensate can be approximated as: 



((A0(r)) 2 ) « Q, 



qua 



r 2 . 



(23) 



Different from the result given by Eq. (21), the phase 
fluctuations are proportional to r 2 when r is much 



smaller than the correlation time r n 



Therefore, if 



Tphasc calculated from Eq. ( [20] ) is much smaller than T qua , 
the analytical result of the phase diffusion time 7~ phaso is 
then 



it fr 



phase 



SN qua dn((N ) ,T) jd (No)' 



(24) 



where SN a 



We now turn to discuss the phase correlation experi- 



ment by the JILA group |21|. The experimental values in 
the experiment are: No = 5 x 10 5 , T « 50 nk, T c « 500 
nk, and a s ~ 5 x 10 -7 cm. It may be helpful to make a 
comparison between the particle number fluctuations due 
to quantum fluctuations and thermal fluctuations. For 
temperature much lower than the critical temperature, 
the analytical result (S 2 N th ) = ir 2 N (T/T c ) 3 /6£ (3) ||] 
can give a rather well description for the particle number 
fluctuations due to thermal fluctuations. For the val- 
ues typical for the experiment by the JILA group 0], a 
simple calculation shows that (S 2 N qua ) / (5 2 N t h) = 62.6. 
Therefore, the thermal fluctuations can be safely omitted 
when the phase diffusion process is investigated for the 
experiment by JILA group Using the formula j20| ) 
(or Eq. (^2|)), the numerical result of the phase diffusion 
time Tpha.se is 119 s, which is much larger than the correla- 
tion time r qua . When obtaining this result, we have used 
the exponential form of the correlation function given by 
Eq. (|l6|). In fact, we can obtain T p h aS e directly from Eqs. 
(0) and (|l4|), and it is worth pointing out that there is 
no important correction to the phase diffusion time, in 
comparison with the result obtained by using the expo- 
nential form (|l6|) . The merit of the exponential form ( |l6| ) 
is that it clearly shows the role of particle number fluctu- 
ations on the phase diffusion process, and the analytical 



result of the phase fluctuations is rather concise using 
this exponential form. 

If the correlation time of the collective excitations is as- 
sumed to be infinity, however, using Eqs. (|l9|) and (g), 
the numerical result of T phaso is 0.62 s, which is much 
smaller than the result given by Eq. (p2|). Although the 
phase fluctuations due to quantum fluctuations are inves- 
tigated in deep in Ref. |2FJ, the finiteness of the correla- 
tion time of the quantum fluctuations was not considered, 
and the dephasing time was approximated as 1 s. In ad- 
dition, it is worth pointing out that although a Langevin 
equation was proposed by Graham |l4],|l5| to investigate 
the phase fluctuations due to thermal fluctuations and 
quantum fluctuations, the phase fluctuations due to col- 
lective excitations were proportional to r 2 , because the 
finiteness of the longevity of the collective excitations was 
not considered too. 

For temperature close to zero, our result of the phase 
diffusion time given by Eq. (^2j) is reasonable because 
of two reasons: (i) In the experiment by the JILA group 
|ilf , the phase of the condensate was found to be very ro- 
bust. In fact, the rigidity of the phase was also shown in 
other experiments, such as the observation of the interfer- 
ence between two BECs ]32||, a nd the recent experiments 
where the optical lat t ic e |33|j3~^ ] is used to investigate 
the coherent properties of the BECs. For example, re- 
cently a BEC Q is created with up to 2 x 10 5 atoms 
and no discernible thermal component. The radial trap- 
ping frequencies are relaxed over a period of 500 ms to 
24 Hz such that the harmonic potential becomes spheri- 
cally symmetric. Then three optical standing waves are 
aligned orthogonal to each other, in order to form a 
three-dimensional lattice potential. In this situation, the 
condensate is distributed over more than 150, 000 lattice 
sites. When the magnetic trap and lattice potential are 
both switched off, it is interested to find that there is 
a high-contrast interference pattern, which means that 
phase is still robust after the BEC has been formed for 
nearly 1 s, and even after the interference between a large 
number of BECs. (ii) In the present work, the correla- 
tion time of the quantum fluctuations is calculated and 
found to be much smaller than the time scale of the phase 
diffusion time. In this situation, we should regard the 
quantum fluctuations as a white noise to investigate the 
phase diffusion process. Recall that the particular collec- 
tive excitations are rather stable when it is created in the 
experiment |29]-|3l|| by applying a small time-dependent 
perturbation, it seems the longevity of the collective exci- 
tations is very long. However, we should note that in the 
problem discussed here for the mechanism of the phase 
diffusion, the collective excitations are created and an- 
nihilated through the quantum fluctuations. As pointed 
out in this paper, the longevity of these collective excita- 
tions is found to be much smaller than the phase diffusion 
time. 

We should note that when obtaining the phase diffu- 
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sion time, the time-energy uncertainty relation is used 
to calculate approximately the longevity of the collec- 
tive modes created due to quantum fluctuations. A more 
accurate average longevity of these collective modes can 
be obtained when interparticle interaction effect is in- 
cluded. Nevertheless, a reasonable order of magnitude 
on the phase diffusion time can be obtained, using Eq. 
( pp| ) and the time-energy uncertainty relation. Addition- 
ally, the collective time r qua can be obtained from Eq. 
( pp| ) when r p h a so is measured in experiment. This gives 
us a chance to check whether standard many-body the- 
ory can be used successfully to investigate the interested 
quantity r qua . 



IV. DISCUSSION AND CONCLUSION 

In summary, the phase diffusion time of the condensate 
due to quantum fluctuations is discussed at extremely low 
temperature. Because the correlation time of the quan- 
tum fluctuations is much smaller than the time scale of 
the phase diffusion, the quantum fluctuations can be re- 
garded as white noise when the phase diffusion of the 
condensate is investigated for extremely low tempera- 
ture. In this situation, the phase diffusion time calcu- 
lated here is much longer than the result obtained in 
the previous theoretical researches [jl4|,|l5],^0| . It is ob- 
vious that the present work can not be applied directly 
to the experiment conducted by the JILA group 
where the complicated rearrangement of the two conden- 
sates would be very important to the phase diffusion pro- 
cess. We shall extend our idea in a future work to the 
case of the specific situation realized by the JILA group 
|l| and discuss the phase diffusion at finite tempera- 
ture. Recently, BECs have been realized in quasi-one 
and quasi- two dimensions [p5| , where new phenomena 
such as quasicondensates with a fluctuating phase [36|-|3Sll 
may be observed. A simple method is developed recently 
p5| , [39| -|4l]| to discuss the particle number fluctuations of 
the low-dimensional condensate. This makes it possible 
to discuss the phase diffusion process in low-dimensional 
condensates. 
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